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ABSTRACT 


This thesis represents the most recent step in the development of a new 
wholeslale level repairable item inventory model for the Navy. Distinguishing 
features of the process by which repairable items are managed in the Navy are 
batching of both repair and procurement orders and a specific repair assessment 
policy which sees repairable items enter repair agencies in batches but leave as 
single items. Approximate steady state probability distributions for inventory 
position and net inventory are developed and these are extensively tested against 
simulation results. Graphical analysis and testing using nonparametric methods 


suggest that the approximations developed closely resemble the distributions 


obtained via simulation. 
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EXECUTIVE SUMMARY 


A. PURPOSE 

The purpose of this thesis is to develop approximately 
correct probability distributions which describe the steady 
State performance of a Navy repairable item inventory. 
Specifically, the distributions of inventory position and net 
inventory are of interest. Once developed these distributions 
can be used in an inventory management mathematical model for 


repairable items. 


B. BACKGROUND 

Since late 1986 work has been undertaken at the Naval 
Postgraduate School to develop a new mathematical model for 
wholesale level repairable item inventory control. The first 
step in the development process is to determine probability 
distributions which describe the steady state behavior of a 
repairable item inventory. The complex mix of stochastic 
processes which occur in the life of a repairable item has 
discouraged previous analytical approaches. This has led to 


the application of simulation techniques to determine 


approximately correct distributions. 





C. RESULTS 

This thesis has shown that it is possible to incorporate 
an analytical approach in the development of approximate 
probability distributions for inventory position and net 
inventory. Using established stochastic modeling results and 
the previously derived approximation for inventory position, 
approximations for the mean and variance of net inventory were 
developed. Two approximating distributions for net inventory 
were proposed and their performance along with that of the 
inventory position approximation were tested against results 
obtained through simulation. 

The results of this thesis support previous research 
results that the convolution of two discrete Uniform random 
variables provides a robust approximation for the steady state 
distribution of inventory position. 

The approximations developed in this thesis for the mean 
of net inventory closely matched simulation output, whereas 
the variance estimates consistently overestimated the 
simulation results. Despite these varied outcomes, the Normal 
and convolution approximations proposed for the steady state 
distribution of net inventory appeared to closely resemble the 
simulation output. On the basis of the graphical analysis 
and testing using nonparametric methods the convolution 
approximation seems a robust approximation for net inventory 


over a wide range of system inputs. 


D. RECOMMENDATIONS 

This thesis represents a significant step in the 
development of a repairable item inventory control model for 
the Navy. The next logical step is for the approximations 
developed in this thesis to be compared against observed 
inventory performance. Subject to a favorable outcome to this 
comparison a readiness-based repairable item inventory model 
using the approximations developed in this thesis should be 


readily forthcoming. 


1. 


I. INTRODUCTION 


A. BACKGROUND 

Significant changes in wholesale supply management 
responsibility have been implemented throughout the Department 
of Defense and the Services since the 1960’s. One of the most 
recent changes has come with the Defense Management Review 
Decision 926, "Consumable Item Transfer", which included the 
transferring of certain stock management responsibilities from 
the Navy to the Defense Logistics Agency (DLA). One effect 
of such changes has been to alter the make-up of wholesale 
level secondary support item [Ref. 1] inventories maintained 
by the Navy. 

Secondary support items used in the Navy are broadly 
classified as either consumable or repairable. In general, 
consumable items are consumed in use and disposed of at the 
time of failure. On the other hand, repairable items can be 
made to function by a repair process after they fail. 

In the past the Navy has managed significant wholesale 
level inventories of both consumable and repairable items. 
However, the changes in supply management responsibilities 
mentioned above have resulted in most consumable items being 
transferred to DLA. Thus, the Navy’s primary focus for 


wholesale level inventory management is now repairable items. 


The differing nature of repairable and consumable items 
dictates that there will be differences in the way that the 
respective inventories are managed. The Navy has separate 
mathematical models which aim to determine optimal strategies 
for managing wholesale level inventories of consumable and 
repairable items. However, the current repairable item model, 
which was developed initially in the late 1960’s and modified 
in 1984, is not supported by any rigorous mathematical 
development and remains founded upon a minor variation of the 
model used to manage consumable items. It is not surprising 
then that the repairable item model displays a low level of 
efficacy. 

Since late 1986 work has been undertaken at the Naval 
Postgraduate School to develop a new mathematical model for 
wholesale level inventory management of repairable items. 
Early attempts to develop an analytical model were 
unsuccessful because of the complex mix of stochastic 
processes which occur in the life of a repairable item. Asa 
result, the focus of recent work has shifted to using 
simulation techniques to study the process. 

The objective of such efforts is to determine 
approximately correct theoretical probability distributions 
which describe the steady state behavior of a repairable item 
inventory. These distributions then can be used in an 


inventory management mathematical model for repairable items. 


B. PRIOR ANALYSIS 

The thesis completed by Maher [Ref. 2] is the most recent 
step in developing an inventory control model for Navy 
repairable items using simulation techniques. While previous 
simulations had been developed, Maher’s work remains 
significant because: 

¢ the simulation model incorporated refinements that more 
accurately describe how wholesale level repairable item 
inventories are maintained in the Navy; 

* some progress was made in the verification of a previously 
derived probability distribution for inventory position, 
and in the identification of a probability distribution 
for net inventory; 

* in the specific context of the Navy repairable item 
inventory, a definition of safety stock was proposed. 

In reviewing the general inventory model literature with 
respect to this thesis four system descriptors are 
significant: 

¢ Batch Ordering policies - which prescribe the quantity of 
a repairable item to be procured or repaired. 

« Repair Assessment/Inspection policies - which prescribe 
the process by which repairable items are initially 
assessed and prepared for repair at a repair center. 

« Shipment Policies - which prescribe the process by which 
repaired and procured items are shipped from a supplier to 


a customer. 


« Condemnation - being the process by which failed items are 
assessed as not being capable of repair. 


As will be described in this thesis, the Navy has specific 
requirements for each of these system descriptors in the 
management of repairable items. These requirements introduce 
complications which are not addressed by the current Navy 
repairable item model or by other prominent repairable item 
inventory models such as METRIC [Ref. 3]. This is not to say 
that there have not been inventory models proposed which deal 
with these more difficult issues. Indeed, the literature has 
examples of batch ordering and/or repair assessment models 
[Ref. 4,5,6, 7 and 8] which, in the main, apply approximating 
techniques to determine optimal inventory control strategies. 
However, the detailed nature of the repairables process does 
not lend itself to the adaptation of such models to the 
specific requirements of the Navy. For example, none of the 
models in the literature consider both repair and procurement 
order batching, or repair induction processes which see failed 
carcasses entering repair agencies as batches but leaving as 
Single items. As a consequence, the Navy’s repairables 
process has been subject to individual and detailed analysis 


at the Naval Postgraduate School. 


Cc. OBJECTIVES 
This thesis follows directly on from Maher’s work with the 
specific objectives to: 
e use simulation output to study the probability 
distributions of net inventory and inventory position for 
a Navy wholesale level repairable item inventory, 
¢ determine the functional relationship between inputs to 
the repairable item management system and parameters of 


these distributions, 


« develop approximately correct theoretical probability 
distributions for net inventory and inventory position. 


D. SCOPE 

Maher simulated a Navy repairable item inventory for a 
specific set of system inputs and employed regression analysis 
techniques to conduct preliminary investigations into 
functional relationships between these inputs and 
distributional parameters. In contrast, this thesis conducts 
the same simulation but over a greater input range, and uses 
stochastic modeling techniques to derive approximations to 


distributional parameters. 


E. PREVIEW 

Chapter II provides a plain language system description of 
the processes which work together to determine the wholesale 
levels of Navy repairable items at any point in time. Chapters 
III and IV develop stochastic models for repairable item 


wholesale inventory levels to support approximations for the 


steady state distribution of inventory position and net 
inventory. Chapter V contains the comparison and analysis of 
results obtained via simulation and stochastic modeling. 
Finally, Chapter VI summarizes the thesis research and 
findings, offers conclusions, and makes recommendations for 


further analysis. 


II. CONCEPTUAL MODEL 


A plain language system description is a key ingredient of 
any study which seeks to understand the relationships between 
various system components and/or £0 predict system 
performance. In the case of this thesis, such a description 
forms the basis for both the simulation of the repairable 
system and the analysis of its stochastic features. Hereafter 
this description shall be referred to as the "conceptual 
model". 

The conceptual model for the management of wholesale level 
inventories of repairable items in the Navy assumes that there 
is a maximum inventory level, SW, and actual inventory levels 
vary from this maximum as a result of demands being placed by 
customers and two distinct but interrelated replenishment 


processes: 


* the repair process 


« the procurement process 


To outline the conceptual model the movement of a 
repairable item between a ship (representing the ultimate user 
of a repairable item), inventory manager, wholesale stock, 


repair agency, and supplier will be described. A summary of 


the system inputs and variables used in the description and 
analysis of the conceptual model is contained in Appendix A. 
The conceptual model assumes that demands for repairable 
items occur at a rate of one unit per requisition and that 
demand follows a Poisson distribution with a mean rate of D 
units per quarter. . On receipt of a demand for a repairable 
item the inventory manager will satisfy the demand from 
wholesale stock on a "first come first served" basis or place 
the demand on backorder. The replenishment of wholesale 
stock levels following these demands will then be undertaken 


by the repair and procurement processes. 


A. THE REPAIR PROCESS 

When a repairable item fails onboard, the ship demands a 
replacement item and dispatches the failed item back to the 
inventory manager so that repair can be arranged. The 
beginning of the repair process for an individual repairable 
item is marked by the return of a failed repairable item to 
the inventory manager. The conceptual model assumes that 
when this return happens it is completed instantaneously. 
However, it is possible that the demand for a replacement item 
ils not accompanied by the successful return of a failed 
repairable item to the inventory manager (for example, due to 


loss of the item in transit between the ship and the inventory 


manager). Therefore, the probability that any one demand will 
be accompanied by a failed item, or carcass, is called the 
Carcass Return Rate (CRR). 

The inventory manager accumulates carcasses until a set 
number (Q,) is reached in the repair batching process. At 
this time Q, failed items are dispatched as a single batch to 
any one of an infinite number of repair agencies. Each 
repair agency is assumed to have infinite capacity to 
undertake repair. Not all carcasses in the batch will be 
capable of repair and thus initially the repair agency 
sequentially assesses and, if necessary, prepares the items 
for repair. The actual decision as to whether a carcass is 
repairable or not takes zero time, whereas the preparation for 
repair, described below, takes a fixed length of time. The 
probability that any one carcass randomly selected will be 
capable of repair is called the Repair Survival Rate (RSR). 
The first carcass is assessed at the moment the batch reaches 
the repair agency. At time t, when a carcass is assessed for 
repair, one of two outcomes will occur. 

If the carcass is assessed as being capable of repair then 
a repaired item will be returned to the wholesale stock at 
time t+RTAT, where RTAT is the known and fixed time interval 
referred to as the Repair Turn Around Time (RTAT). The repair 
assessment of the next item in the batch (should there be one) 
will then occur at time t+REP, where REP is a known and fixed 


time interval. If the carcass is assessed as not being 


capable of repair then the inventory manager is immediately 
notified, the quantity in repair reduced by one, and the 
"attrition" is recorded in the procurement process at time t. 
The repair assessment of the next item in the batch (should 
there be one) will then occur instantaneously at time t. 
Therefore the time element REP can be thought of as the repair 
preparation time or the time delay in commencing actual repair 
caused by the need to conduct diagnostic tests, fault 
isolation procedures and material preparation etc for an 
individual repairable item prior to considering the next item 
in the batch. On the basis of the repair process described 


above REP is necessarily less than or equal to RTAT. 


B. THE PROCUREMENT PROCESS 

When a ship demands a replacement for a failed repairable 
item but, for whatever reason, the failed item is not returned 
to the inventory manager, or an item is assessed in the repair 
process as not being capable of repair then an attrition is 
said to have occurred. Attritions require the procurement of 
additional items to maintain inventory levels and, thus, an 
attrition marks the start of the procurement process. As 
with the repair process, the inventory manager accumulates 
attritions until a set number (Q,) is reached in the 
procurement batching process. At this time the inventory 
manager places an order of size Q, with a supplier which 


arrives in wholesale stock (as a single batch) a fixed time 


10 


interval later. This time interval is referred to as the 
Procurement Lead Time (PCLT) and is assumed to be 
significantly greater than RTAT. 

Figure 1 shows the movement of repaired and procured items 
(RFI - Ready for Issue), carcasses and attritions in the 


repair and procurement processes. 





Orders at 
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—_ RFI Items 
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== Attritions 
Repalr 


Batching 










Batching 





Figure 1. Item Flow Chart for Repairable Item Conceptual 
Model. 


C. MAGNITUDE OF MODEL INPUTS 

Based on the conceptual model described, there are eight 
system inputs which will determine the distribution of 
inventory levels at any particular point in time. 

The range of system inputs that have been considered in 
this thesis together with the basis for their selection are 
outlined as follows: 

* Quarterly Demand Rate (D): In the case of repairable items 
the quarterly demand rate tends to be low. For example, 


on the basis of sample information provided by Mr. John 
Boyarski at the Ships Parts Control Center, demand rates 


less than 9 units per quarter tend to dominate. Rates 
extending to 40 units per quarter are, however, 
considered. 


¢ Carcass Return Rate (CRR): Financial redemption values 
applied to repairable items in the Navy provide incentives 
for users to return failed carcasses to the inventory 
manager. Return rates in the range from 0.4 to 1.0 are 


considered. 

¢ Repair Survival Rate (RSR): Repairing a failed item 
consumes considerable effort and resources. As a 
consequence, this would dictate that failed carcasses have 
a reasonable probability of being repaired. Repair 
survival rates in the range from 0.4 to 1.0 are 
considered. 


¢ REP: At the component and sub-assembly levels it seems 
reasonable to expect that this time element might be 
reasonably short, perhaps not exceeding two weeks. In 
more complex sub-assembly and assembly levels this delay 
due to diagnostic testing and fault isolation, etc., may 
be significant. REP delays in the range from 0.0 to 1.6 
quarters are considered. 


¢ Repair Turn Around Time (RTAT): Historical evidence shows 
that RTAT varies between organic and commercial repair 
agencies. On the basis of data provided by Mr. John 
Boyarski turnaround times in the range from 0.3 to 2.3 
quarters are considered. 
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¢ Procurement Lead Time (PCLT): Repairable items generally 
experience long resupply times. On the basis of 
information provided by Mr. Boyarski lead times in the 
range from 3.5 to 8.5 quarters are considered. 


* Batch/Order Quantities (Q,/0,): While classified as inputs 
for the purposes of this thesis, in any inventory model 
these quantities will become decision variables. 
Therefore, in testing distributional approximations these 
inputs should remain reasonably unconstrained. However, 
there are some practical relationships that might be 
reasonably expected among Q,, D, REP and RTAT. Repairable 
items that experience low demand rates and/or large 
RTAT/REP values are likely to be batched for repair in 
small quantities. The converse is also true. Therefore 
Batch/Order quantities up to 42 units which apply these 
relationships are considered. 


* Maximum Inventory Level (SW): The nature of the conceptual 
model is such that any value of SW would serve equally 
well for testing purposes. To provide continuity the 
same value as used by Maher, 72 units, is used throughout 
this thesis. 

The stochastic modeling conducted in this thesis also 
assumes that there is some level of material flow through the 
procurement and repair processes. This assumption translates 
to the joint condition that CRR>0 and product CRR*RSR<1. The 
case when CRR equals zero represents the strict procurement 
model, the Navy application of which was discussed by Maher 
[Ref. 2:pp. 6-18]. The case when CRR*RSR equals one 
represents the strict repair model. While this later 
situation is highly unlikely to arise, the arguments developed 


in this thesis can be easily extended to cover such a case. 


1.3 


D. MEASURES OF INTEREST 
The inventory level measures which are of primary interest 
to this thesis are Inventory Position (IP) and Net Inventory 
(NI) and are defined as follows: 
IP = On Hand RFI Stock + On Order + In Repair - 
Backorders, 


which, in the context of the conceptual model, is equivalent 


to 
IP = SW - Carcasses and attritions currently being 
accumulated by the Inventory Manager. 
NI = On Hand RFI Stock - Backorders, 


which is equivalent to, 
NI = IP - On Order - In Repair. 
Here RFI refers to the items which are "as new" and Ready for 


Issue. 
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III. INSTANTANEOUS REPAIR ASSESSMENT 


The case of instantaneous repair assessment is when the 
time element REP in the conceptual model is considered 
insignificant and set to zero. Under this condition, as a 
repair batch of size Q, reaches a repair agency, the entire 
batch will be assessed instantaneously. Thus, all carcasses 
commence repair immediately and attritions are processed 
immediately. This case simplifies to some degree the 
stochastic analysis of the conceptual model, while at the same 
time providing a foundation for the more difficult scenario of 
non-instantaneous repair assessment. 

The first section of this chapter reviews the development 
of the steady state distribution of inventory position. The 
second section considers the net inventory level for the case 
of instantaneous repair assessment, with key results being an 
approximation for the mean and variance, leading to a 


suggested probability distribution for net inventory. 


A. STEADY STATE DISTRIBUTION OF INVENTORY POSITION 

Maher [Ref. 2:pp. 26-32] outlined the development of the 
steady state distribution of inventory position for the 
repairables model as previously derived by Professor McMasters 


at the Naval Postgraduate School. To remain consistent with 


LS 


the notation introduced in the conceptual model, and for the 
purposes of clarity, the supporting arguments and result shall 
be summarized in this thesis. 

To introduce notation, let 

I(t) = the inventory position at time t; 

BR (t) = the. number of carcasses in the repair 


batching process at time t; and 


ii] 


B,(t) the number of attritions in the procurement 
batching process at time t. 


Then, it follows from the conceptual model that 
I(t) = SW- B,(t)- B,(t). (1) 


As demands are generated and accompanied by carcasses, 
Ba(t) will increase by one until Q, items are held in the 
repair batching process. At this later time a repair order 
will be generated instantaneously decreasing the number in the 
repair batching process to zero. It follows from the 
conceptual model that this demand will occur as a Poisson 
process, and therefore the repair batching process can be 
modeled as a continuous time Markov chain with state 
space{0,1,2 ... Q.-1} and a doubly stochastic transition 
Matrix [Ref 9:p. 191]. Thus, the steady state distribution 


for the number of carcasses batching is given by 


P(B,(t) =w) = for w=0,1,2,...-Qp-1. (2) 


R 
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Similarly, as attritions are generated, B,(t) will 
increase by one until Q, items are held in the procurement 
batching process. At this later time a procurement order 
will be generated instantaneously decreasing the number in the 
procurement batching process to zero. As stated in the 
conceptual model, these attritions result both from demands 
not being accompanied by a carcass and carcasses being 
assessed beyond repair in the repair process. It follows 
from the conceptual model that this first source of attritions 
will form a Poisson process and occur at a constant rate. 
Results from Maher’s work support the assumption that the 
second source of attrition can be approximated as a Poisson 
process [Ref 2.pp. 29-30]. Thus, the procurement batching 
process can also be modeled as a continuous time Markov chain 
with state space {0,1,2 ...Q,-1} and a doubly stochastic 
transition matrix. The steady state distribution for the 
number of attritions batching is approximated by 

P(B,(t)=v) # Z for v = 0,1,2,.+.Q,-1. (3) 

There is no intuitive reason why at a particular point in 
time the two batching processes would be dependent and, 
therefore, the steady state probability distribution for 
inventory position can be approximated from the convolution of 
P(Ba(t)) and P(Be(t)). The resulting distribution for P(I(t)) 


is shown in equation (4). 
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xt1 





for 0<x<x, 


0,0 
Min(Qpz, Op) 
P[I(t)= SW-x] = |————" "2 for x,<x< 
O:0p Se (4) 

Xyaxt1-X 

— for <XSX, 

0;0z Or Xx, ‘Max 
0 otherwise; 


where Xy,x,, X:, and x, are defined as 


Xyax = OptQp-23 
xX, = Min(Q,,Qp) -1; (5) 
% * Xyax~X,: 


Given that inventory position is approximated utilizing the 
convolution of two discrete Uniform random variables, it 


follows that 


Qp-1 a Qp-1 . 


= 6 
E[I(t)] = sw s a7 (6) 





zs (Q,-1) (Qp+1) + (Q,-1) (Q,+1) 


Var[Z(t)] 3 2 


(7) 


B. STEADY STATE DISTRIBUTION OF NET INVENTORY 
To consider the distribution of net inventory let 


N(t) = net inventory level at time t; 


D(t,,t>) number of demands in the period (t,,t,); 


D’ (t,,t2) 


number of demands in the period (t,,t,) 


accompanied by a carcass; 
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DB (hyve = number of demands in the period (t,,t>) 

not accompanied by a carcass; and 

Reg lg te) = number of carcasses inducted into repair 

agencies in the period (t,,t,) that will 
be assessed as repairable. 

Since the time any item spends at a repair agency (which 
is at most RTAT) is strictly less than PCLT, it follows that 
in the time period (t-PCLT,t) 

¢ Any quantity in procurement or actual repair at time t- 
PCLT will have been received by time t. 

« Any repairable carcasses commencing actual repair in the 
period (t~PCLT,t-RTAT) will have been received by time t. 
With instantaneous repair assessment, for any repairable 

carcass to commence actual repair in the time period (t- 
PCLT, t-RTAT) it must have been inducted into a repair agency 


in the time period (t-PCLT,t-RTAT). Thus it follows that 
N(t) = I(t-PCLT) - D(t-PCLT,t) + Ry,(t-PCLT, t-RTAT), (8) 


which, after substituting equation (1) for I(t-PCLT), and 
noting that the Poisson random variable D(t,,t,) can be 


decomposed into the sum of D’(t,,t,) and D"(t,,t,), becomes 


N(t) = SW - Bp(t-PCLT) + Ry(t-PCLT,t-RTAT) - D!(t-PCLT, t) 
- B,(t-PCLT )- D!(t-PCLT,t). 
(9) 
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1. Repairable Items (R,,) Inducted Into Repair Agencies 
As carcasses enter the repair process they are 
initially batched until a quantity Q, is reached, at which 
point a repair order is instantaneously placed. From the 
conceptual model it is known that the inter-arrival times of 
these carcasses into the repair process are independent and 
identically distributed Exponential random variables with 
parameter (D*CRR). Thus, the times between the dispatching of 
successive repair batches of size Q, to a repair agency will 
be independent and identically distributed Erlang random 
variables with parameters (Q,,D*CRR). Let 
Tos = the time between the dispatching of successive 
repair batches of size Q, to a repair agency; 


and so 





Og 
= A (10) 
E( Typ] D CRR’ 
Var [T,_] = ——2 _. (11) 
(D CRR)? 
In addition, let 
O,(t) = the number of repair batches of size Q, inducted 


into repair agencies in the time interval [0,t). 
The counting process {0,(t), t>0} is a Renewal process as a 
consequence of the independent and identically distributed 
times between dispatching of successive repair batches. 
In addition, from elementary Renewal theory the mean 


value function (i.e., the expected number of repair batches of 
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size Q, sent to repair agencies in the time interval (t,,t,)) 


can be approximated using Blackwell’s theorem [Ref 10:p.41] as 


(6-6) D CRR (t,-t,) 
BIO 4E;) Ey] 3 = ee 
FS ae he BLT, —— eo 


the number of repair batches of size Q, 
inducted into repair agencies in the time 
period (t,,t,). 
Furthermore, if 
Wi = the number of repairable carcasses in the ith 
repair order of size Qa, 
then w,; are independent and identically distributed Binomial 


random variables with mean and variance defined by 


Elw,] = Op RSR; 


Var[w,] = Op RSR (1-RSR). 


Therefore the counting process {R,,(0,t), t>0} 
classified as a Renewal Reward Process [Ref 10:p. or 


Cumulative Process [Ref 1l:p. 91] described by 


Or( Cy, C2) 
Wy £Or ONG, ty) = Vy 250% 


Rip(t,,0,) = > 


(¢) otherwise. 





The approximate expected number of items inducted into 
repair agencies in the period (t,,t,) that will be assessed 


as repairable follows as 


EIR, (t,t)! 


Bl ela, tt, #,) fo,ft,, £31 ] 
E[O,(t,,t,) Q, RSR] 


(16) 


Op RSR E[O,(t,, t,)] 


1] 


D CRR RSR (t,-t,). 


Conditioning can also be applied to approximate the 
variance of R,,(t,,t,). Relaxing the time interval notation, 


then 


Var[R,p] = E[Var[R,,/0,]] + Var[E[R,,/0,]] 
E[O, Op RSR (1-RSR)] + Var[0O, Op RSR] (17) 


Qp RSR (1-RSR) E[O,] + Q,? RSR? Var[O,]; 


where the variance for the Renewal Process {0,;(t), t>0} is 


approximated by [Ref 11:p.70] 


Var [Top] (t,-t,) 
E[T ppl? 

, DCRR (t,-t;) 

OR 


Var[O,(t,,t,)] * 
(18) 


Substitution of equations (12) and (18) into (17) leads to 


Var [Rip(t,,t,)] = D CRR RSR (t,-t,). (19) 
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A further result of interest is the covariance between 
O, and Ryg, which follows from Cox [Ref 11:p.94] as 


Cov [Rrzp,O,] = Var[O,] Op RSR 


_ DCRR RSR (t,-t,) (20) 
Or , 


2. Expected Net Inventory 
The preceding results can be used to determine the 
expected net inventory. Specifically, taking the expected 


value of equation (9) and relaxing the time interval notation, 


E(N(t)] = SW - E[Bp] + E[Ry] - E[D’] 








(21) 
- E[B,] - E[D”). 
Using the results of equations (6) and (16) leads to 
=i 
E([N(t)] = Sw - oe" + D CRR RSR(-RTAT+PCLT) 
pit (22) 
- D CRR PCLT - ar - D(1-CRR) PCLT; 
which simplifies to 
win) «en- 24. SS 
2 2 
(23) 


- D CRR RSR RTAT 
- D(1-CRR RSR) PCLT. 


By way of validation, consider the situation where Q, 
and Q, are both set equal to one, and REP is set equal to 
zero. In this situation a carcass arriving at the repair 


process at time t is immediately sent to a repair agency to 
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undergo instantaneous assessment and, should it be an 
attrition, it will enter the procurement process at the same 
time t. Thus, with infinite repair and supply capacity this 
situation is the same as operating two independent infinite 
server queues with Poisson arrivals and fixed service times. 
Specifically, this situation will be characterized by 
¢ An infinite server queue in the repair process 
experiencing Poisson demand arrivals with rate D*CRR*RSR 
and fixed service time of RTAT; 
¢e An infinite server queue in the procurement process 
experiencing Poisson demand arrivals with rate 
D*(1-CRR*RSR) and fixed service time of PCLT. 


From Queuing Theory [Ref 9:pp. 225-226] it is known that the 


expected net inventory level in this case will be 


E([N(t)] = SW - D CRR RSR RTAT - D (1 - CRR RSR) PCLT. (24) 


The same result follows from substitution of Q, and OQ, equal 


to one in equation (23). 
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3. Variance of Net Inventory 
To consider the random variable dependencies in the 
case of instantaneous repair assessment a convenient way to 


rewrite equation (9) is 


N(t) = SW- R(t) - R’(t) - P(t); (25) 
where 
R(t) = Bp(t-PCLT) - R,,(t-PCLT,t~-RTAT) + D!(t-PCLT, t-RTAT) ; 
R!(t) = D/(t-RTAT, t); 
P(t) = B)(t-PCLT) + D"(t-PCLT,t). 


(26) 

The independence of the repair batching process (B,) 
and the procurement batching processes (B,) underlies the 
development of the steady state distribution of inventory 
position. In addition, established results for the Poisson 
process dictate that D" and D’ are independent, being Poisson 
processes formed from the categorization of the Poisson 
process D. Finally, R,,(t-PCLT,t-RTAT) and B,(t-PCLT) are 
formed from demand type D’ which occurs before time t-RTAT and 
t-PCLT, respectively. Thus, being functions of independent 
random variables, it follows that R(t), R’(t) and P(t) are 
independent. Furthermore, the number of attritions in the 


procurement batching process at a specific point in time is 
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independent of the demand which follows in a _ subsequent 


period. Therefore, B,(t-PCLT) and D"(t-PCLT,t) are 
independent. Based on these facts 

Var(N(t)] = Var[R(t)] + Var[R/(t)] + Var[P(t)]; (27) 
where 


w= [@p ~ 1) (Op * 1), a rR) Pcrr; 
= (28) 


Var[R’/(t)] = D CRR RTAT. 


Var [P(t)1 


To examine the variance of R(t), consider the repair 
batching process at time t-PCLT. This process will continue 
to accumulate carcasses until a total of Q, is held. In this 
case the time elapsed subsequent to t-PCLT will be equivalent 
to the length of time from t-PCLT until the next renewal in 
the Renewal Process formed by the generation of repair orders. 
In the situation when this remaining batching time is large in 
comparison to the time period (t-PCLT,t-RTAT), any repairable 
item inducted into a repair agency in the period (t-PCLT,t- 
RTAT) is highly likely to have been in the repair batching 
process at t-PCLT. Thus, under these circumstances a 
Significant degree of dependency would be expected between 
Ba(t~PCLT) and R,,(t-PCLT,t-RTAT) . An alternative situation 
is when the remaining batching time is small in comparison to 
PCLT-RTAT. In this case any repairable item inducted into a 
repair agency in the time period (t-PCLT,t-RTAT) is more 


likely to be as a result of demand D’ subsequent to t-PCLT 
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than to have been in the batching process at t-PCLT. These 
two cases illustrate the expected dependencies between B,(t- 
PCLT), R,g(t-PCLT,t-RTAT) and D’ (t-PCLT,t) under instantaneous 
repair assessment. Additionally, based on these arguments the 
covariance between B,(t-PCLT) and R,,g(t-PCLT,t-RTAT) is 
expected to be non-negative. 

Thus, relaxing time interval notation, it follows that 


the variance of R(t) is 


Var[R] = Var[B,] + Var[R,,] + Var[D’] 
(29) 
~ 2 Cov[By, Ry] - 2 Cov[Ryp,D’). 


To determine the covariance between R,;, and D’ consider that 
in any fixed time interval (t-PCLT,t-RTAT) 
D!(t-PCLT, t-RTAT) = Q, O,(t-PCLT, t-RTAT) + a; (30) 


where @ is a random variable in the range 


0x a@ 5 Q,-1 for O,=0; 
(31) 
-(Qp-1) < @ < (Q,-1) for O,21. 
Thus 
Cov[Ryp,D!] = Cov[Ryp, 0,0,+@] . (32) 


In a probabilistic sense, the number of repairable carcasses 
entering repair agencies, R,,, is a function of the number of 
repair batches of size Q, accumulated by the Inventory 
Manager, O,;, and the Repair Survival Rate, RSR. Therefore it 


seems reasonable that R,;,(t-PCLT,t-RTAT) is independent of @. 
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Combining equations (30) and (32) with equation (20), together 


with the independence assertion, leads to 


Cov[Ryp,D’] = Cov[Rzp,0,0;] + Cov[Rzp, &] 


O, Cov[Rzp, O;] (33) 


n 


D CRR RSR(PCLT-RTAT) . 


Therefore, under conditions where the covariance between B, 
and R,, is insignificant, the substitution of equations (19) 
and (33) into equation (29) leads to 


= (p71) (Opt) _ _ 
5 DCRR RSR(PCLT-RTAT) (34) 


+ D CRR(PCLT-RTAT) ; 


Var [R(t)] 


and the variance of net inventory equation (27) becomes 


3 (Op-1) (O,+1) 4 (Op-1) (O,+1) 
12 12 (35) 


+ D CRR RSR RTAT + D(1-CRR RSR) PCLT. 


Var [N(t)] 


Once again from elementary Queuing Theory results, the 
variance for the boundary case where Q, and Q, are both set 
equal to one and REP is set equal to zero (being the operation 
of two independent infinite server queues with Poisson 
arrivals and fixed service times) is provided by equation 
(35): 

4. Probability Distribution for Net Inventory 

The task now is to identify a probability distribution 

which characterizes the steady state levels of net inventory 


in the case of instantaneous repair assessment. One way to 
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identify such a distribution is to match moments of the 
approximating distribution with those determined through 
stochastic analysis of the conceptual model. For lower order 
moments, a necesSary condition must be that the mean and 
variance of the chosen distribution are reasonably well 
approximated by equations (23) and (35), respectively. 

After conducting a number of simulation runs for a 
specific range of conceptual model inputs, Maher [Ref 2:pp. 
65-67] found strong evidence to suggest that the Normal 
distribution was a good approximating distribution for net 
inventory. Thus, one approach is to model N(t) as a Normal 
random variable with mean and variance given by equations (23) 
and (35), respectively. A Normal approximation is likely to 
allow algebraic simplification in the subsequent development 
of a readiness based inventory model for repairables. 
However, it also carries with it the condition that the 
probability that net inventory is larger than SW is strictly 
positive; a condition that violates the conceptual model. 

A second approach to approximating the distribution of 


net inventory is to let net inventory be described by 
N(t) = SW- [Bg(t) + By(t) + ¥(t)], (36) 


where Ba(t) and B,(t) are as described previously, and Y(t) is 


a Poisson random variable with mean wp defined as 


B = DCRR RSR RTAT + D(1-CRR RSR) PCLT. (37) 


Zo 


Then the approximating distribution for N(t), based on 
equation (36), would be the convolution of three independent 
random variables B,(t), Bp(t) and y(t), and would have its mean 
and variance described by equations (23) and (35), 
respectively. In this approach, since Ba(t), B,(t) and Y(t) 
are all positive random variables, the probability that N(t) 
ls greater than SW is strictly zero. 

The probability density function for this convolution 
is 


1 we (j+1) pee oe eo 
20, & S51 or 0<ZSx, 





1 x (j+1) p7Je# 
0,0p > (z-J) ! 


ap min (Op, Op) = p2Je “# 


OO, yin (2-9)! cai eae 





P[N(t) = SW-z] = 


1 wh (j+1) p2%e 
0,0p j=0 (z-j) ! 





, min (Op, Qa) =a pzJe# 
0,0, jouer (2-9)! 


2 Cea 7) pe 


L 0,0, JRX2+1 (z-7J) ! 





for X,<Z; 


(38) 


where Xy,,, X,, and x, were defined by equation (5). 
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Iv. NON-INSTANTANEOUS REPAIR ASSESSMENT 


This chapter considers the case where the time element REP 
in the conceptual model assumes some positive value. Under 
this condition, as a repair batch of size Q, reaches a repair 
agency, each item in the batch is sequentially assessed before 
proceeding into repair or being declared an attrition. Thus, 
carcasses enter repair agencies in batches but they are 
assessed individually and, if repaired, returned to wholesale 
stock as individual items. 

Generally speaking, repairable item inventory models 
discussed in the literature assume that this repair assessment 
happens instantaneously [Ref. 6 and 8] (as presented in the 
previous chapter), and/or that repair orders sent to repair 
agencies are assessed as a single batch and returned as a 
batch after repair [Ref 5,6,7 and 8]. Such conditions combined 
with manageable batch arrival distributions obviously have 
resulted in tractable results. However, the method of repair 
assessment and return of repaired units outlined in the 
conceptual model and described above, together with the 
batching of repair orders initially incorporated in Chapter 
III, distinguish the Navy model from these other repairable 


models discussed in the literature. 
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The nature of the conceptual model (with the condition of 
non-instantaneous repair assessment) has not made it possible 
in this study to develop a strong analytical argument leading 
to the probability distribution of net inventory. However, 
building on the arguments presented for the instantaneous 
repair assessment case, an approximation for the mean of net 
inventory can be devised and a simplistic approach applied to 
approximate the distribution. 

The first section of this chapter briefly considers the 
steady state distribution of inventory position. The second 
section considers the net inventory level for the case of non- 
instantaneous repair assessment, leading to a suggested 


probability distribution for net inventory. 


A. INVENTORY POSITION 

Maher [Ref 2:p.59] found that for the range of system 
parameters he considered the variation of REP had little 
effect on the simulated inventory position distribution. Lt 
remained essentially that stated by equation (4). Thus, a 
starting point in this analysis is to choose to model 
inventory position in the same way as was done in the case of 
instantaneous repair assessment. Specifically, the steady 
state distribution of inventory position in the case of non- 


instantaneous assessment is assumed to be approximated by 


32 


equations (4) and (5). The results of this thesis to be 
presented later support Maher’s results and this general 


approach for modeling inventory position. 


B. NET INVENTORY LEVELS 
To consider net inventory levels let N(t), D(t,,t,), 


D’(t,,t,), D"(t,,t.), and R,,(t,,t,) be as defined previously, 


and let 
A, (t) = number of attritions previously inducted into 
repair agencies and still awaiting repair 
assessment at time t; and 
A, (t) = number of repairable carcasses previously 


inducted into repair agencies and still 
awaiting repair assessment at time t. 

As with the case of instantaneous repair assessment, the 
general conditions apply that under the assumption that RTAT 
is strictly less than PCLT in the time period (t-PCLT,t) 

* Any quantity in procurement or actual repair at time t- 

PCLT will have been received by time t. 

« Any repairable carcasses commencing actual repair in the 

period (t-PCLT,t-RTAT) will have been received by time t. 
Accounting for the fact that carcasses entering repair 


agencies now may be delayed by the repair assessment action 
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before being declared repairable or attritions, it follows 


that 


N(t) = I(t-PCLT) - D(t-PCLT,t) - A,(t-PCLT) 
(39) 
+ Rip(t-PCLT, t-RTAT) - A,(t-RTAT) ; 


which, after substituting equation (1) for I(t-PCLT) becomes 


N(t) = SW - Bp(t-PCLT) + Ryp(t-PCLT, t-RTAT) - D(t-PCLT, t) 
- A,(t-PCLT) - A,(t-RTAT) - B,(t-PCLT ). (40) 


It should be noted that when REP is set to zero then A, and A, 
are necessarily zero, and equation (40) reduces to equation 


(9) « 


1. Repairable Carcasses (A,) and Attritions (A,) At 
Repair Agencies Awaiting Repair Assessment 


The robust results of Little [Ref. 12] can be used to 
study the number of repairable carcasses and attritions at 
repair agencies awaiting repair assessment. Little showed 


that in a steady-state queuing process the following formula 


holds 
L= Aw; (41) 
where 
L = expected number of units in the system; 


= expected time spent by a unit in the system; and 


= expected rate of arrivals to the system. 
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The robustness of this result lies in the observation that it 
is 
remarkably free of specific assumptions about arrival 
and service distributions, independence of inter-arrival 
times, number of channels, queue disciplines, etc. [Ref. 
12:p.387] 
It was stated previously that the number of repairable 
carcasses entering repair agencies by time t {R,g(0,t), t>0} 
can be modeled as a Reward Renewal Process or Cumulative 


Process. In particular, Renewal Theory results for this type 


of process lead to 
E(R,p(0,t)] = DCRR RSR t, (42) 


as derived in equation (16) in Chapter III. Thus, choosing 
to approximate the steady state rate of repairable carcasses 
entering repair agencies as D*CRR*RSR and defining 
De = time a randomly selected item spends at a repair 
agency awaiting repair assessment, 


it follows from Little’s formula that 


E[Az(t)] = DCRR RSR E[T,]. (43) 


Similar arguments can be made for the number of 


attritions entering repair agencies by time t such that 


E[A,(t)] = DCRR (1-RSR) E[T,]. (44) 
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2. Time Awaiting Repair Assessment (T,) 
To determine the time any carcass or attrition spends 
at a repair agency awaiting repair assessment (T,), define the 


random variables 


ne = the sequence number of an item in a repair batch 
of size QO, (1; 2... Q,)s; and 
Re. JS the number of items in a particular repair batch 


of size Q, assessed as repairable prior to an item 
undergoing repair assessment (0, 1... Q,-1). 

The conceptual model states that as a repair batch of 
size Q, enters a repair agency the first item in the batch is 
assessed immediately. The outcome of this assessment 
directly influences the amount of time the remainder of that 
batch spends waiting for repair assessment to commence. In 
contrast, the assessment of the last item in the batch does 
not influence the time any item spends awaiting repair 
assessment. More generally, the time any randomly selected 
item spends waiting for repair assessment to commence will be 
an integer sum of time packages REP and will depend on the 
number of repairable items in the batch of size Q, that are 


assessed before it. This can be represented as the product 


T, = R, (REP). (45) 
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Then, conditioning on the sequence number of an item, the 


conditional probability becomes 
P(R,=x/I=i) = ( #7) (rsp) *(1-RsR)4* for 0sxsi-1; (46) 
where any position in the order is equally likely; i.e., 
P(iei) = 4, (47) 
Op 


Then the expected time any item spends at a repair agency 


awaiting repair assessment will be determined as follows: 





E[T,] = REP E[R,]; (48) 
E(R,] = El ELR,/I) 1; (49) 
E[R,/I=i] =(i-1) (RSR) ; (50) 
EIR iy (4-1) RSR P(I=i) = RsR £271, 
A i (51) 
=1 
E[T,]=£2— rsp Rep. (52) 





2 
3. Expected Net Inventory 
Relaxing the time interval notation and taking the 


expectation of N(t) given by equation (40) leads to 


E([N(t)] = Sw - E[B,] + E[R,,] - E[D] 
- E[A,] - E[A,] - E[B,). 


(53) 


37 


Substituting the results obtained from equations (6), (16), 


(43) and (44) leads to 


-1 -1 
Oet _ 2h _ a _oRR RSR) PCLT 
2 2 (54) 


- D CRR(RSR RTAT + E[T,]). 


E([N(t)] = Sw - 


Finally, adding the result given by equation (52) gives 


Qp-1 QOp-1 








E(N(t)] = Sw ~ ~2— - “2 - D(1-CRR RSR) PCLT 
(55) 
- D CRR RSR(RTAT + —— ; 


It should be noted that when REP is set to zero equation (55) 
reduces to equation (23) in Chapter III. 
4. Variance of Net Inventory 

The structure of N(t) given by equation (35) allows 
for the direct application of the arguments used in the 
previous chapter to the determination of variance of net 
inventory in the case of non-instantaneous repair assessment. 

However, point estimates for the variances of A, and A,, and 
the form of any dependencies induced by these random 
variables, are not obvious. 

While formulae for higher order moments similar to 
Little’s result do exist in the case of Exponential inter- 
arrival patterns [Ref. 13], they are ineffective in this case 
as in the conceptual model, Erlang inter-arrival patterns of 
repair batches at an infinite number of repair agencies are 


experienced. Then, on arrival at a repair agency, each batch 
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enters a Single server queue for repair assessment. The 
inability to accurately model this process leads to 
consideration of the following naive analysis. 

In the boundary cases where the Chapter III assumption 
concerning the insignificance of Cov[B,,R,;,] 1s applied, Q, is 
set equal to one, and/or REP is set equal to zero, any 
expression for the variance of net inventory should equate to 
equation (35). Similarly, as REP approaches zero the point 
estimate results should approach those for the instantaneous 
repair assessment case. Furthermore, for REP sufficiently 
small, it might be expected that the lower order moments for 
the deviation of net inventory below SW approach those of a 
random variable formed from the convolution of Ba(t), Bp(t) and 


a Poisson random variable with mean 


= D (1-CRR RSR) PCLT 


(0,71) REP 


+ D CRR RSR(RTAT + 3 


ie 


With this naive analysis the variance of net inventory would 


be estimated by 


ro (Op-1) (Op+1) % (Qp-1) (Op+1) 
12 12 


Var[N(t)] 


(O,-1) coil 


+ DCRR RSR (RTAT + 3 


+ D(1-CRR RSR) PCLT. 





5. Probability Distribution for Net Inventory 

Once again using the approach of matching lower order 
moments, and applying the results of the naive analysis, an 
approximation can be made for the steady state distribution of 
net inventory. 

Maher found. that the Normal distribution served as a 
good approximation for net inventory for both instantaneous 
and non-instantaneous repair assessment. Thus, one approach 
1s to model N(t) as a Normal random variable with mean and 
variance given by equations (55) and (57), respectively. 

A second approach, identical to that followed in 


Chapter III, is to let net inventory be described by 


N(t) = SW - [BR(t) + Bp(t) + ¥(t)]; (58) 


where B,(t) and B,(t) are as describe previously, and Y(t) is 
a Poisson random variable with mean given by equation (56). 
Then the approximating distribution for N(t) resulting from 
equation (58) and the convolution of three independent random 
variables Ba(t), B,(t) and Y(t) will have mean and variance 
described by equations (55) and (57), respectively. The 
probability distribution function for this convolution will be 


that given by equation (38) with wf now given by equation (56). 
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Vv. RESULTS 


The purpose of this chapter is to discuss the performance 
of the approximations developed for inventory position and net 
inventory in Chapters III and IV. The first section of the 
chapter outlines the methods used to test the distributional 
approximations developed previously. The second section of 
the chapter discusses results obtained in the instantaneous 


repair assessment case. Finally, the third section discusses 


the results for the case of non-instantaneous repair 


assessment. 


A. TESTING METHOD 

The objective of this thesis is to develop probability 
distributions which describe the steady state behavior of the 
Navy's wholesale repairable item inventory. These 
distributions can then be used in an inventory management 
model for repairable items. Typically such models seek to 
optimize a selected measure of effectiveness such as total 
variable cost of maintaining inventory, supply response time, 
material availability, etc. Thus, the most appropriate form 
of testing would be to investigate how well the suggested 
probability distributions work in such a model. In the 


absence of an inventory management model a more stringent test 





is to investigate how well the approximate analytical 
probability distributions for inventory position and net 
inventory compare with empirical distributions generated 
through simulation. If the results are satisfactory then the 
distributions can be used in the development of an inventory 
management model. This later approach was adopted in this 
thesis. 

To validate the probability distributions for inventory 
position and net inventory which were described in Chapters 
III and IV, a number of simulation runs were conducted using 
the MODSIM II high-level programming language. MODSIM II is 
a general-purpose, modular, object-oriented programming 
language which has built-in capability to perform discrete 
event simulation [Ref 14: pp.1-2]. The MODSIM II simulation 
implements the conceptual model described in Chapter II as was 
that previously studied by Maher [Ref. 2] using SIGMA 
(SImulation Graphical Modeling and Analysis system). The 
simulation output published by Maher therefore provided an 
initial verification tool for the MODSIM II program used in 
this thesis. 

For the range of system inputs considered, and program run 
times experienced, an appropriate simulation length was 
considered to be 1000 quarters. Based on analysis of 
simulation time series an initial transient time of 50 


quarters was generally applied. 
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The simulation runs were used to study the performance of 
point estimates and distribution approximations using a 
combination of graphical techniques and goodness of fit tests. 
As with Maher’s thesis, for the Chi-square goodness of fit 
test the "observed frequencies" resulting from the simulation 
model were determined by comparing the total length of time 
that each value of the inventory position and net inventory 
occurred in steady state with the total run time less initial 
transient time (i.e., the time-weighted distributions were 
determined). As Maher noted [Ref. 2:p. 58], this use of the 
Chi-square goodness of fit test stretches the theory and is 
not exact. However, no other goodness of fit test exists for 


this type of frequency determination. 


B. INSTANTANEOUS REPAIR ASSESSMENT 
1. Point Estimates 

As an initial illustration of the effectiveness of the 
approximations developed for the particular case of 
instantaneous repair assessment a number of system input 
settings were considered. These settings fall within the 
range of system input values discussed in Chapter II, and are 
shown in Table 1. 

Point estimates for the mean and variance of inventory 
position and net inventory for these five system input 
settings are shown in Tables 2 and 3, respectively. [In all 


cases it can be seen that the approximations are very close to 
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those obtained from simulation runs. In particular, over a 
large number of additional runs the simulated results gave a 
mean and variance of inventory position and a mean of net 
inventory which were always very close to the approximate 
formulae. Simulation output for the variance of net 


inventory tended to-be less close to the approximate formulae. 





Table 1 System Input Settings 





2. Steady State Distribution of Inventory Position 
A graphic comparison of the simulated and 
approximating probability distributions for inventory position 
resulting from the application of input settings in Table 1 is 
contained in Appendix B. These initial results support 
Maher’s findings [Ref 2:p.78] that the approximation developed 


provides a very good fit. 
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Table 2. Comparison of Inventory Position Point Estimates. 





Table 3. Comparison of Net Inventory Point Estimates. 


Appendix C contains results of Chi-Square goodness of 
fit test for 78 simulation runs applying a range of system 


inputs. In all cases the observed test statistic, x’, is less 
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than the critical value, x’... Thus, at the 5% significance 
level the conclusion is to fail to reject the null hypothesis 
that the distribution of inventory position is as described by 
equations (4) and (5). 

3. Steady State Distribution of Net Inventory 

Graphical comparisons of the simulated distributions 
for net inventory and Normal approximations resulting from the 
application of input settings in Table 1 to the simulation 
model are shown in Appendix D, and the - convolution 
approximations are shown in Appendix E. In these initial 
cases it is difficult to discern any significant difference 
between the graphics for the Normal and convolution 
approximations. Both seem to provide a very good 
representation of the simulated output. However, it should 
be noted that the graphics for the Normal approximation do not 
include theoretical probabilities for net inventory levels 
greater than SW (72 in all cases) because net inventory levels 
will never exceed SW in the conceptual model. 

Appendices F and G contain the Chi-Square goodness of 
fit test results for 78 simulation runs using the Normal and 
convolution approximations, respectively. Once again at the 5% 
level of significance the general conclusion is to fail to 
reject the null hypothesis that the distribution of net 
inventory is described by the Normal approximation with mean 


and variance given by (23) and (35), respectively, or by the 
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convolution given by equation (38). Tt is only in the cases 
where low variability of net inventory is experienced in 
steady state (reflected in the degrees of freedom for the Chi- 
square test which in all cases are equivalent to the number of 
inventory levels experienced in steady state less one) that 
the Normal approximation appears to consistently fail. This 
failure is likely to be the result of the Normal approximation 
applying significant positive probability to net inventory 
levels greater than SW. The convolution approximation is 
accepted in all cases. 

The variance approximation derived in Chapter III and 
applied in cases 1 through 5 above has assumed that the 
covariance between B,(t-PCLT) and R,,(t-PCLT,t-RTAT) is 
insignificant. Given the relationship developed in equation 
(29), and the intuitive reasoning developed in Chapter III 
supporting positive correlation between these two random 
variables, it follows that the variance approximation should 
be greater than the variance reflected in simulation results. 
Indeed, this phenomena was consistently observed over the 
large number of simulation runs performed during this study. 
Despite this result, the convolution approximation for the 
steady state distribution of net inventory seems to work very 
well and is supported by the Chi-square goodness of fit test 


results. 
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C. NON-INSTANTANEOUS REPAIR ASSESSMENT 
1. Steady State Distribution of Inventory Position 

Appendix H contains the Chi-square goodness of fit 
results for the approximating distribution of net inventory 
following the application of a range of system input settings 
to the simulation model. In these cases REP varies from 0.0 
to: 0.5. Once again the observed test statistic yx’ is 
continually, and significantly, less than the critical value, 
X°.05.0- Thus, at the 5% significance level the conclusion is 
to fail to reject the null hypothesis that the distribution of 
inventory position is as described by equations (4) and (5). 

These results together with those obtained in the 
instantaneous repair assessment case add further weight to 
Maher’s findings; that the convolution of two discrete Uniform 
random variables provides a robust approximation for the 
steady state distribution of inventory position. 

2. Net Inventory Point Estimates 

To assess the approximations developed for net 
inventory the theoretical point estimates were first compared 
to simulation output using the input settings in Table 4. 
Graphical comparisons of the theoretical and simulated mean 
and variance of net inventory are at Appendices I and J, 
respectively. Each input setting was replicated four times 


in the simulation model. 
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Table 4. System Input Settings. 


As with the case of instantaneous repair assessment, 
the theoretical estimates for the expected net inventory are 
very close to those gained from the simulation. However, in 
the case of the variance the theoretical point estimate 
consistently overstates the result obtained through 
simulation. This result is not surprising given the 
simplistic basis for the variance estimate. The simulation 
results also show that variance as a distributional parameter 
is not particularly well behaved. The nature of this behavior 
presents difficulties in trying to determine a general 
estimate for the difference between the simulated and 


estimated variances using regression techniques. 
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3. Steady State Distribution of Net Inventory 

Despite the shortcomings displayed in the variance 
estimates for net inventory, the Normal and convolution 
distributional approximations were tested against the 
simulation model. Appendices K and L contain the Chi-Square 
goodness of fit test results for 90 simulation runs applying 
a range of system input settings and using the Normal and 
convolution approximations, respectively. Once again at the 
5% level of significance the general conclusion is to fail to 
reject the null hypothesis that the distribution of net 
inventory is described by the convolution approximation given 
by equation (38). The Normal approximation, however, 
continues to consistently fail the hypothesis test where low 
variability of net inventory is experienced. Note that, as 
expected from the conceptual model, in the case where Q, 
equals one the value of REP has no influence on the 
distributional parameters. 

The simulation runs conducted in Appendices K and L 
with REP equal to 0.5 (its largest value in Appendices K and 
L) were used to provide a graphical comparison of the net 
inventory probabilities obtained through simulation and the 
convolution approximation as given by equation (38). These 
plots are contained in Appendix M, with case input settings 
summarized in Table 5. The plots show that the convolution 


resembles the simulated result quite closely, with the 


50 


absolute differences between the simulated and approximated 
probabilities being very small in all cases. While the plots 
for Cases 18 and 19 appear to show that the two distributions 
begin to differ, it should be noted that in these cases the 
relative shape of the distribution function is extremely flat 
and the absolute differences in probabilities remain very 
small. The Chi-square test for these two cases (see last two 
results in Appendix L) leads to the conclusion, at the 5% 
Significance level, that the convolution is a reasonable 
representation of the steady state distribution of net 


inventory. 





Table 5. System Input Settings. 
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vi. SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 


A. SUMMARY 

In the past the Navy has managed substantial inventories 
of both consumable and repairable secondary items. However, 
changes implemented in the Department of Defense wholesale 
supply management aS a consequence of Defense Management 
Review Decision 926 have seen a large migration of the 
consumable items to DLA and leaves the Navy and the other 
services to focus their inventory management efforts on 
repairable items. In addition, readiness-based sparing is 
becoming increasingly important as a result of Congressional 
pressures. To assist in this adjustment of focus the Navy 
requires tools which help in the determination of optimal 
strategies for the maintenance and control of a repairable 
item inventory. 

The Navy’s current repairable item inventory control model 
has been found to give unsatisfactory results when attempts 
are made to use it for readiness-based sparing. As a 
consequence, there is a requirement for a fundamental analysis 
of the Navy’s repairable item inventory management process, 


rather than a simple patch on the existing model. In light of 
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this, work has been undertaken at the Naval Postgraduate 
School to develop a new mathematical model for wholesale level 
repairable item inventory control. 

The first step in the development process is to determine 
probability distributions which describe the steady state 
behavior of a repairable item inventory. The complex mix of 
stochastic processes which occur in the life of a repairable 
item has discouraged previous analytical approaches. This 
has led to the application of simulation techniques to 
determine approximately correct distributions. 

This thesis has shown that it is possible to use the 
simulation results to suggest and to test an analytical 
approach for the development of approximate probability 
distributions for inventory position and net inventory. The 
repairable item inventory process was analyzed for the two 
separate cases of instantaneous and non-instantaneous repair 
assessment. Using established stochastic modeling results and 
the previously derived approximation for inventory position, 
approximations for the mean and variance of net inventory were 
developed. Two approximating distributions for net inventory 
were proposed and their performance along with that of the 
inventory position approximation were tested against results 


obtained through simulation. 


53 


B. CONCLUSIONS 
1. Inventory Position 

The results of this thesis support Maher’s findings 
[Ref.2] that the convolution of two discrete Uniform random 
variables provides a robust approximation for the steady state 
distribution of inventory position. This result holds for 
both the instantaneous and the non-instantaneous repair 
assessment cases over a wide range of system inputs. 

2. Net Inventory 

The approximations developed in this thesis for the 
mean of net inventory closely matched simulation output in 
both the instantaneous and non-instantaneous repair assessment 
cases. The variance estimates consistently overestimated the 
Simulation results, however. A major reason for this appears 
to be the assuming away of some random variable dependencies 
(i.e., covariances were assumed to be zero). Despite these 
results, the Normal and convolution approximations proposed 
for the steady state distribution of net inventory appeared to 
closely resemble the simulation output. On the basis of the 
graphical analyses and goodness of fit testing the convolution 
approximation seems a robust approximation for the net 


inventory distribution over a wide range of system inputs. 
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C. RECOMMENDATIONS 

This thesis represents a significant step in the 
development of a repairable item inventory control model for 
the Navy. The next logical step is for the approximations 
developed in this thesis to be compared against inventory 
characteristics observed for actual repairable items being 
managed by the Navy. In particular, a formula for the 
expected number of time-weighted backorders at any instant of 
time needs to be developed and the values of average days 
delay it computes compared to those being observed by the 
Navy’s Inventory Control Point repairable item inventory 
managers. Subject to a favorable outcome to this comparison 
a readiness-based sparing repairable item inventory model 
using the approximations developed in this thesis should be 


readily forthcoming. 
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APPENDIX A. SUMMARY OF SYSTEM INPUTS AND RANDOM VARIABLES 


System Inputs 


SW 


RSR 


Qe 
Qe 


REP 


RTAT 


PCLT 


Maximum inventory level 
Quarterly demand rate 


(Carcass Return Rate) Probability that a demand will 
be accompanied by a carcass 


(Repair Survival Rate) Probability that a carcass 
(failed item) is repairable 


Number of items in a repair batch/order 
Number of items in a procurement batch/order 


The time interval that elapses in the repair 
assessment process following the assessment of a 
repairable item until the assessment of the next 
carcass in the order 


(Repair Turn Around Time) The time taken by a repair 
agency to repair an individual item including repair 
assessment time 


(Procurement Lead Time) The time taken by a supplier 
to fill a procurement order 


Random Variables 


Tt) 


Ba (t) 


Bp (t) 


N(t) 


Ditty, ty) 


I 


the inventory position at time t 


= the number of carcasses in the repair 
batching process at time t 


= the number of attritions in the procurement 
batching process at time t 


= net inventory level at time t 


= number of demands in the period (t,,t,) 
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DY (tt, , 6b) 
I (ty 7 Ge) 
Rip (t,, to) 
Tre 

O,(t) 
O,(t,,t2) 


number of demands in the period (t,,t,) 
accompanied by a carcass 


number of demands in the period (t,,t,) not 
accompanied by a carcass 


number of carcasses inducted into repair 
agencies in the period (t,,t,) that will be 
assessed as repairable 


the time between the despatching of 
successive repair batches of size Q, to a 
repair agency 


the number of repair batches of size Q, 
inducted into repair agencies in the time 
interval [0,t) 


the number of repair batches of size Q, 
inducted into repair agencies in the time 
period (t,t) 


Random Variables - Non-Instantaneous Repair Assessment Only 


Ra 


number of attritions previously inducted into 
repair agencies and still awaiting repair 
assessment at time t 


number of repairable carcasses previously 
inducted into repair agencies and still 
awaiting repair assessment at time t 


time an item spends at a repair agency 
awaiting repair assessment 


the sequence number of an item in a repair 
batch of size Q, (1, 2 ... Qg) 


the number of items in a particular repair 
batch of size Q, assessed as repairable prior 
to an item undergoing repair assessment (0, 1 
«+s Qp-d) 
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APPENDIX B. COMPARISON OF PROBABILITY DISTRIBUTIONS FOR 
INVENTORY POSITION WITH INSTANTANEOUS REPAIR ASSESSMENT 


INSTANTANEOUS ASSESSMENT 
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Probability 


INSTANTANEOUS ASSESSMENT 
Case 2 
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INSTANTANEOUS ASSESSMENT 


Case 3 
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Probability 


INSTANTANEOUS ASSESSMENT 
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APPENDIX C. GOODNESS OF FIT TEST RESULTS FOR INVENTORY 
POSITION WITH INSTANTANEOUS REPAIR ASSESSMENT 


D CRR RSR  RTAT PCLT Q, Q x? ) ree v 
1 0.56 0.86 0.9 8.21 1 1 0 0 0 
1.2 0.93 0.97 2.03 3.8 2 2 4.376772 5.99147 2 
1.4 0.57 0.88 1.23 7.48 3 3 4.210447 9.48773 4 
1.6 0.4 0.86 1.21 5.7 3 4 0.875907 1.0705 5 
1.8 0.46 0.91 0.98 4.07 4 7 7.337342 16.919 9 
2 0.7 0.62 0.7 4.55 3 2 0.82123 7.81473 3 
2.2 0:55. 4 223 8.04 4 4 2.489711 12.5916 6 
2.4 0.55 0.94 0.99 7.43 3 6 4.671086 14.0671 7 
2.6 0.67 0.65 0.6 3.99 2 6 3.48185 12.5916 6 
2.8 0.81 0.73 2.08 7.83 4 8 8.605134 18.307 10 
3 0.72 0.94 0.98 4.05 4 7 6.688924 16.919 9 
3.2 0.77 0.68 1.79 7.61 4 10 3.268578 21.0261 12 
3.4 0.54 0.48 2.26 5.04 5 10 3.55471 22.3621 13 
3.6 0.74 0.68 1.51 6.61 7 8 4.930263 22.3621 13 
3.8 0.77 0.93 1.21 4.28 6 9 3.899164 22.3621 13 
4 0.64 0.42 1.17 5.04 8 8 5.562638 23.6848 14 
4.2 0.83 0.44 0.54 6.11 4 6 3.731272 15.5073 8 
4.4 0.65 0.94 2.27 4.57 7 7 5.118714 21.0261 12 
4.6 0.91 0.9 qa Pie’ 6 15 7.146105 30.1435 19 
4.8 0.9 0.98 0.38 7.94 5 20 9.674624 35.1725 23 
5 0.67 0.51 0.31 5.25 5 12 4.146568 24.9958 15 
5.2 0.93 0.6 2.06 7.71 7 24 12.79699 42.5569 29 
5.4 0.76 0.6 0.51 4.61 8 10 3.918378 26.2962 16 
5.6 0.7 0.5 2.24 5.57 6 6 4.317124 18.307 10 
5.8 0.82 0.98 0.53 7.09 7 6 2.271654 19.6751 11 


6 0.73 0.94 1.58 6.76 8 14 6.151449 31.4104 20 





So oo 2 ©: 
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Qe 


10 


12 


10 


10 


11 


13 
14 


10 


20 


10 
17 
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x? 
3.82707 
1.821077 
6.411615 
1.49633 
4.506006 
12.42674 
1, 782152 
3.042209 
2.121086 
4.087285 
3.196558 
1.778286 
3.109613 
1.351535 
1.522632 
3.074357 
0.772842 
2.276205 
3.276276 
2.632506 
2.23805 
2.456715 
1.217982 
2.163939 
2.144016 
23956159 
2.204993 
4.160763 
1.078626 
0.696505 
2.047428 


18. 


46. 


0.97 


38 


34 


65 


2 


x 


-58592 

- 832667 
- 866691 
-576058 
-724438 
-900257 
+ 962103 
-118684 
-285073 
-504164 
-767047 
-079228 
«275526 
-503216 
- 008606 
-548498 
- 919398 
-B12295 
-660321 
-296568 


-627907 


APPENDIX D. COMPARISON OF PROBABILITY DISTRIBUTIONS FOR 
NET INVENTORY WITH INSTANTANEOUS REPAIR ASSESSMENT - NORMAL 
APPROXIMATION 
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Probability 


INSTANTANEOUS ASSESSMENT 


0.16 





Case 3 - Normal Approximation 
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Probability 


INSTANTANEOUS ASSESSMENT 


Case 4 - Normal Approximation 
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Probability 


INSTANTANEOUS ASSESSMENT 


Case 5 - Normal Approximation 
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APPENDIX E. COMPARISON OF PROBABILITY DISTRIBUTIONS FOR 
REPAIR ASSESSMENT = 


NET INVENTORY WITH INSTANTANEOUS 
CONVOLUTION APPROXIMATION 


INSTANTANEOUS ASSESSMENT 


Case 1 - Convolution Approximation 


Probability 
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INSTANTANEOUS ASSESSMENT 


Case 5 - Convolution Approximation 








Net Inventory 
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APPENDIX F. GOODNESS OF FIT TEST RESULTS FOR NET INVENTORY 
WITH INSTANTANEOUS REPAIR ASSESSMENT - NORMAL APPROXIMATION 


D CRR RSR RTAT PCLT Q, Q> 0 ce v 
1 0.56 0.86 0.9 8.21 1 1 36.437592 22.3621 13 
1.2 0.93 0.97 2.03 3.8 2 2 32.577991 16.919 9 
1.4 0.57 0.88 1.23 7.18 3 3 34.632537 24.9958 15 
1.6 0.4 0.86 1.21 5.7 3 4 29.968999 31.4104 20 
1.8 0.46 0.91 0.98 4.07 4 7 20.7532 30.1435 19 

0.7 0.62 0.7 4.55 3 2 20.661384 28.8693 18 
2.2 0.55 4 23 8.04 4 4 23.959783 40.1133 27 
2.4 0.55 0.94 6.99 7.43 3 6 36.810676 32.6705 21 
2.6 0.67 0.65 0.6 3.99 2 6 20.667706 28.8693 18 
2.8 0.81 0.73 2.08 7.83 4 8 16.28019 40.1133 27 
3 0.72 0.94 0.98 4.05 4 7 11.198242 33.9244 22 
3.2 0.77 0.68 41.79 7.61 4 10 28.682596 44.9854 31 
3.4 0.54 0.48 2.26 5.04 5 10 18.828359 43.7729 30 
3.6 0.74 0,68 1.517 6.61. 7 8 46.356123 46.1942 32 
3.8 0.77 0.93 1.21 4.28 6 9 12.333993 41.3372 28 
4 0.64 0.42 1.17 5.04 8 8 23.19932 43.7729 30 
4.2 0.83 0.44 0.54 6.11 4 6 10.994421 42.5569 29 
4.4 0.65 0.94 2.27 4.57 7 7 14.800585 49.9854 31 
4.6 0.91 0.9 Gif 7.5 6 15  27.655204 44.9854 31 
4.8 0.9 0.98 0.38 7.94 5 20 58.173473 47.3998 33 
5 0.67 0.51 0.31 5.25 5 12 13.423188 49.8019 35 
5.2 0.93 0.6 2.06 7.71 7 24 38.561637 74.4686 56 
5.4 0.76 0.6 0.51 4.61 8 10 19.506471 52.1925 37 
5.6 0.7 0.5 2.24). F572 16 6 15.481147 48.6023 34 
5.8 0.82 0.98 0.53 7.09 7 6 8.338465 43.7729 30 
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APPENDIX G. GOODNESS OF FIT TEST RESULTS FOR NET INVENTORY 


WITH INSTANTANEOUS REPAIR ASSESSMENT - CONVOLUTION 
APPROXIMATION 

D  CRR  RSR- RTAT PCLT Q, Q, x V ss. v 
1 0.56 0.86 0.9 -8.21 1 1: 12.761395 22.3621 13 
1.2 0.93 0.97 2.03 3.8 2 2 14.084955 16.919 9 
1.4 0.57 0.88 1.23 7.18 3 3. 18.074717 24.9958 15 
1.6 0.4 0.86 1.21 5.7 3 4 -.21.024362 31.4104 20 
1.8 0.46 0.91 0.98 4.07 4 7  11.444031 30.1435 19 
2 0.7 0.62 0.7 4.55 3 2 8.504929 28.8693 18 
2.2 0.55 1 2.3 8.04 4 4 12.622673 40.1133 27 
2.4 0.55 0.94 0.99 7.43 3 6 25.873481 32.6705 21 
2.6 0.67 0.65 0.6 3.99 2 6 14.631932 28.8693 18 
2.8 0.81 0.73 2.08 7.83 4 8  17.534952 40.1133 27 
3 0.72 0.94 0.98 4.05 4 7 7.619581 33.9244 22 
3.2 0.77 02:68 1.79 F.61 4 10 23.695209 44.9854 31 
3.4 0.54 0.48 2.26 5.04 5 10 16.042827 43.7729 30 
3.6 0.74 0.68 1.51 6.61 7 8 40.4659 46.1942 32 
3.8 0.77 0.93 1.21 4.28 6 9 6.986432 41.3372 28 
4 0.64 0.42 1.17 5.04 8 8  21.608591 43.7729 30 
4.2 0.83 0.44 0.54 6.11 4 6  10.867147 42.5569 29 
4.4 0.65 0.94 2.27 4.57 7 7  17.302609 49.9854 31 
4.6 0.91 0.9 1.11 7.5 6 15 21.098001 44.9854 31 
4.8 0.9 0.98 0.38 7.94 5 20 9.657205 47.3998 33 
5 0.67 0.51 0.31 5.25 5 12 12.253753 49.8019 435 
5.2 0.93 0.6 2.06 7.71 7 24 18.665303 74.4686 56 
5.4 0.76 0.6 0.51 4.61 8 10 19.559415 52.1925 37 
5.6 0.7 0.5 2.24 5.57 6 6 14.471527 48.6023 34 
5.8 0.82 0.98 0.53 7.09 7 6 7.403263 43.7729 30 
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APPENDIX H. GOODNESS OF FIT TEST RESULTS FOR INVENTORY 
POSITION WITH NON-INSTANTANEOUS REPAIR ASSESSMENT 


D CRR RSR REP RTAT PCLT Q, QQ 7 G58, 8 v 
1 0.56 0.86 0.05 0.9 8.21 1 5} 0 ) 0 
2 0.7 0.62 0.05 0.7 4:55 3 2 1.832042 7.81473 3 
3 0.72 0.94 0.05 0.98 4.05 4 7 6.339916 16.919 9 
4 0.64 0.42 0.05 1.17 5.04 8 8 5.942194 23.6848 14 
6 0.73 0.94 0.05 1.58 6.76 8 14. 4.979643 31.4104 20 
8 0.53 0.43 0.05 2.17 6.12 8 20 3.203709 38.8852 26 
12 0.43 0.81 0.05 0.94 3.64 14 9 1.7425 32.6705 21 
16 0.43 0.8 0.05 1.14 7.59 24 20 2.967227 58.1241 42 
20 0.86 0.91 0.05 1.51 5.43 18 22 2.058203 53.3837 38 
1 0.56, 0.86 0.1 0.9 8.21 2 1 0 0 0 
2 0.7 0.62 0.1 0.7 4.55 3 2 0.205652 7.81473 3 
3 0.72 0.94 0.1 0.98 4.05 4 7 5.178484 16.919 9 
4 0.64 0.42 0.1 117 85504 8 8 4.478211 23.6848 14 
6 0.73 0.94 0.1 1.58 6.76 8 14. 2.832181 31.4104 20 
8 0.53 0,43 0.2 BIL? 6.12 8 20 3.506728 38.8852 26 
12. 0.43 «0.81 40.1 0.94 3.64 14 9 2.99248 32.6705 21 
16 0.43 0.8 0.1 1.14 7.59 24 20 3.66927 58.1241 42 
20 0.86 0.91 0.1 1.51 5.43 18 22 1.843208 53.3837 38 
if 0.56 0.86 0.15 0.9 21. J 1 0 0 0) 
2 0.7 0.62 0.15 0.7 4/55 3 2 0.404025 7.81473 3 
3 0.72 0.94 0.15 0.98 4.05 4 7 3.354933 16.919 9 
4 0.64 0.42 0.15 1.17 5.04 8 8 4.789292 23.6848 14 
6 0.73 0.94 O.15 1.58 6.76 6 14 3.338997 31.4104 20 
8 0.53 0.43 0.15 2.17 6.12 8 20 5.382016 38.8852 26 
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APPENDIX I. COMPARISON OF POINT ESTIMATES FOR EXPECTED NET 
INVENTORY WITH NON-INSTANTANEOUS REPAIR ASSESSMENT 
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APPENDIX J. COMPARISON OF POINT ESTIMATES FOR VARIANCE OF 
NET INVENTORY WITH NON-INSTANTANEOUS REPAIR ASSESSMENT 
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APPENDIX K. GOODNESS OF FIT TEST RESULTS FOR NET INVENTORY 


WITH NON- INSTANTANEOUS REPAIR ASSESSMENT = NORMAL 
APPROXIMATION 

D CRR RSR REP RTAT PCLT Q, Q, x? NP oie v 

1 0.56 0.86 0.05 - 0.9 8.21 1 1 36.437592 22.3621 13 
2 Om 0.62 0.05 0.7 4.55 3 2 23.792937 28.8693 18 

3 0.72 0.94 0.05 0.98 4.05 4 7 8.536619 33.9244 22 

4 0.64 0.42 0.05 1.17 5.04 8 8 28.44658 47.3998 33 

6 0.73 6,94 0.05 1,58 6.76 ‘8 14 24.441525 52.1925 37 

8 0.53 0.43 0.05 2.17 6.12 8 20 16.896291 70.9934 53 
12 0.43 0.81 0.05 0.94 3.64 14 9 13.487128 66.3386 49 


16 0.43 0.8 0.05 1.14 7.59 24 20 £42.759145 113.1449 90 


20 0.86 0.91 0.05 1.51 5.43 18 2a 42.027846 75.6238 57 
1 0.56 0.86 0.1 0:9 B.21 1 1 36.437592 22.3621 13 
2 0.7 0.62 O.1 0.7 4.55 38 2 30 .89962 26.2962 16 
3 0.72 0.94 O.1 0.98 4.05 4 7 24.814517 37.6525 25 
4 0.64 0.42 O.1 1.17 5.04 8 8 17.448078 46.1942 32 
6 O73 0.94 O.1 1.58 6.76 8 14 27.785375 52).1925 37 
8 0.53 0.43 O.1 2017 6:12. 8 20 15.899915 70.9934 53 
12 0.43 0.81 O.1 0.94 3.64 14 9 27.586662 65.1706 48 
16 0.43 0.8 0.1 1.14 7.59 24 20 62.235448 90.5312 70 
20 0.86 0.91 O.1 1.51 5.43 18 22 43.938006 76.7774 58 


1 0.56 O.86 0.25 0.9 6.21 1 1 36.437592 22.3621 13 
2 0.7 0.62 0.15 0.7 4.55 3 2 28.671046 26.2962 16 
3 Ore O594 0:15 0.98 4.05 7 14.69815 35.1725 23 
4 0.64 0.42 OFS Jl SeO04 8 8 24.651764 46.1942 32 
6 0.73 0.94 0.15 1.58 G.76 8 14 26.460645 54.5724 39 
8 0.53 O43 O10 2.17 6.12 8 20 21.857007 75.6238 57 
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APPENDIX L. GOODNESS OF FIT TEST RESULTS FOR NET INVENTORY 


WITH NON-INSTANTANEOUS REPAIR ASSESSMENT - CONVOLUTION 
APPROXIMATION 

D CRR RSR REP RTAT PCLT Q, Q be Di ee v 

1 0.56 0.86 0.05 -0.9 G20 4 1 12.761395 22.3621 13 

2 0.7 Os62 0.05 0.7 4.55 3 2 10.580141 28.8693 18 

3 0.72 0.94 0.05 0.98 4.05 4 7 8.990693 33.9244 22 

4 0.64 0.42 0.05 1.17 5.04 8 8 31.889717 47.3998 33 

6 0.73 0.94 0.05 1.58 6.76 8 14 26.063501 52.1925 37 

8 O53 €.43 @Q.05 2.17 6.12 68 20 12.659687 70.9934 53 
12 0.43 0.81 0.05 0.94 3.64 14 9 13.428029 66.3386 49 


16 0.43 0.8 0.05 1.14 7.59 24 20 44.022047 113.1449 90 


20 0.86 0.91 0.05 Lett 5.43 18 ZA 43.118399 75.6238 Ga 
1 0.56 0.86 0.1 0.9 8.21 a: 1 12... 761395 22.3621 18 
Zz Oia? 0.62 0.1 0.7 4.55 3 2 6.87995 26.2962 16 
3 0.72 0.94 0.1 0.98 4.05 4 7 13.431709 37.6525 25 
4 0.64 0.42 0.1 ind? 5.04 8 8 15.213948 46.1942 ae 
6 0:..73 0.94 G.4 1.58 6.76 8 14 26.134285 5201925 37 
8 0:53 0.43 8.1 Aol? 6.12 8 20 14.543465 70.9934 53 
12 0.43 0.81 O52 0.94 3.64 14 9 27.68042 65.1706 48 
16 0.43 0.8 0.1 1.14 7.59 24 20 63;.3203 90.5312 70 
20 0.86 On. FL 0.1 1w5Sl 5.43 18 Zid 45.482399 76.7774 58 
1 0.56 0.86 O.cs 8.9 8.21 1 i 12:. 761395 22:..3621 i3 
2 0.7 0.62 O.15 0.7 4.55 3 2 12.08988 26.2962 16 
3 0.72 0.94 0.15 0.98 4.05 4 7 9.82721 35.1725 Zo 
4 0.64 0.42 O.4de 4.19 5.04 8 8 24.017196 46.1942 32 
6 0.33 0.94 0.05 L.58 6.76 8&8 14 24.261553 54.5724 a3 
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APPENDIX M. COMPARISON OF PROBABILITY DISTRIBUTIONS FOR NET 
INVENTORY WITH NON- INSTANTANEOUS REPAIR ASSESSMENT = 
CONVOLUTION APPROXIMATION 
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